Boron in copper: a perfect misfit in the bulk and cohesion enhancer at a grain 
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Our ab initio study suggests that boron segregation to the E5(310)[001] grain boundary should 
strengthen the boundary up to 1.5 ML coverage (15.24 at/nm^). The maximal effect is observed at 
0.5 ML and corresponds to boron atoms filling exclusively grain boundary interstices. In copper bulk, 
B causes significant distortion both in interstitial and regular lattice sites for which boron atoms are 
either too big or too small. The distortion is compensated to large extent when the interstitial and 
substitutional boron combine together to form a strongly bound dumbell. Our prediction is that 
bound boron impurities should appear in sizable proportion if not dominate in most experimental 
conditions. A large discrepancy between calculated heats of solution and experimental terminal 
solubility of B in Cu is found, indicating either a sound failure of the local density approximation 
or, more likely, strongly overestimated solubility limits in the existing B-Cu phase diagram. 
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Introduction 

Boron has an extremely good record in improving in- 
tergranular cohesion in metals. It is mostly famous for 
curing the long standing problem of room temperature 
brittleness in NisAli Boron segregation was found to 
reinforce grain boundaries in other intermetallic com- 
pounds (FeAl, NiAI, NiaSi) and to improve low temper- 
ature ductility in bcc iron and refractory metals, such as 
Mo and W (see, e.g., Ref. 0] and references therein). 

The effect of boron addition on copper is far less stud- 
ied. According to the Cu-B phase diagram^- boron 
solubility in copper is low, 0.06 at.% at room temper- 
ature rising to 0.29 at.% at the eutectic temperature 
1013 °C. Dissolved boron has a strong propensity to 
segregate to surfaces and interfaces. It is not clear 
whether segregation weakens or strengthen grain bound- 
aries. Nevertheless, doping copper with boron is found 
to be efficient in preventing segregation of antimony to 
grain boundaries.^ Substantial improvement of mechani- 
cal properties of nanocrystalline Cu samples is reported 
as B segregation can be used to limit grain growth during 
heat treatment;^ Despite such encouraging experimental 
findings, quite surprisingly, no theoretical simulations of 
boron at copper grain boundaries or even free surfaces, 
are known to us. 

Similarly poor is the situation with studies of boron 
in bulk copper. It is not even clear whether boron oc- 
cupies interstitial or substitutional positions. Analysis of 
boron's neighbours in the Periodic Table does not rule 
out either possibility. Carbon is an interstitial impurity 
in Cu)^ whereas Al and Be are substitutional impurities. 
In recent experimental work in which accelerated boron 
ions were implanted in Cu, both types of impurities were 
observed 

In the present study we employ standard density func- 
tional calculations to study the behaviour of boron impu- 
rities at a copper grain boundary and in the bulk. We find 
that boron strengthens the E5{310}[001] symmetric tilt 
grain boundary in the whole range of boundary coverages 



investigated (up to 1.5 ML). The maximum strength- 
ening occurs at 0.5 ML at which boron exclusively oc- 
cupies grain boundary interstices. We further identify 
mechanisms responsible for grain boundary strengthen- 
ing within the framework of the "ghost impurity cycle" 
proposed in our previous work.^ The cycle admits the oc- 
cupation of both substitutional and interstitial positions 
by impurity atoms at the interface, a feature that is fully 
exploited in the present study. 

Due to the peculiar interplay of atomic sizes, intersti- 
tial and substitutional positions are equally unwelcoming 
to boron in bulk Cu. Boron is too big an interstitial im- 
purity and too small a substitutional impurity. As a re- 
sult, both have nearly the same heat of solution with the 
interstitial position marginally preferred. If, however, B 
atoms combine in dumbells then most of the elastic dis- 
tortion of the host is eliminated and significant lowering 
of the heat of solution is achieved. The energy gain is 
so large that boron dumbells should persist up to high 
temperatures. Another surprising finding of our study 
is a large disagreement between the theoretical heat of 
solution and experimentally observed maximal solubility 
of B in Cu. 

The paper is organised as follows. SecUoutlines defini- 
tions and thermodynamical relations used in the present 
work. The computational setup is described in Sec. [TTl 
The results on B in bulk Cu are presented in Sec. Illli 
whereas Sec. IIVI is concerned with the effect of boron 
at the grain boundary. In the latter section, we first 
look at the change of the work of separation and at the 
impurity segregation energies (Sec. IIV A[) . then we dis- 
cuss atomic structure of the boundaries and free surfaces 
with different boron content (Sec. IIV 5)) . and finally apply 
the "ghost impurity cycle" to reveal the mechanisms re- 
sponsible for the cohesion enhancement (Sec. IIV C|) . Our 
main findings are summarised in Sec. [V] A thermody- 
namic model used to estimate the concentration of dif- 
ferent kinds of boron impurity in bulk copper, is outlined 
in the Appendix. 
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FIG. 1: "Ghost impurity cycle" used for discussion of the 
effect of impurity on grain boundary strength. A and B de- 
note the pure and segregated grain boundaries, respectively, 
in their equilibrium geometry. Boundary C is created by sub- 
stituting impurity atoms in B (black spheres) with vacancies 
(white spheres) without further atomic relaxation. Configu- 
ration D is A in which the host atoms that will be replaced 
by impurity in B are removed with other atoms kept in place. 
As and Bs denote the free surfaces into which grain bound- 
aries A and B cleave. Cg and Ds are prepared from Bs 
and Cs respectively, using the above strategy, namely, im- 
purity atoms in Cs, or host atoms to be replaced by impu- 
rity in Ds, are removed keeping the positions of other atoms 
fixed. Path A— »D— *C-^B refers to substitutional impuri- 
ties, whereas path A^C^B applies to interstitial impurities 
(Ref.i). 



crease ^cicav Bi in Cu is a classic example of the for- 
mer. Copper grain boundaries with bismuth segregating 
eventually reach the condition ^cicav < ^disi leading to a 
ductile-to-brittlc transition. 

In the limit of fast separation which we assume in 
the present study, any impurity exchange between newly 
formed surfaces and bulk during the decohesion is pre- 
vented. In this limit, t/cicav can be identified with the 
reversible work of separation 

^cleav = Wscp = ^{6^-03"}, (1) 

where A is the surface area, G^^ is the excess Gibbs free 
energy of a representative piece of material containing 
grain boundary, and is the sum of two excess Gibbs 
free energies corresponding to surfaces formed after deco- 
hesion. Eq. ([1]) assumes that the impurity excess F at the 
grain boundary is equal to the sum of surface impurity 
excesses. The excess in Eq. ([T]) is defined with respect to 
the underlying bulk crystal. 

Our ab initio calculations refer to the zero temperature 
limit hence the Gibbs free energies are replaced with total 
energies. It is convenient then to express the changes in 
Wsop due to impurity in terms of segregation energies per 
impurity atom, Egeg- (Segregation energy is the energy 
required to remove all impurity from an interface and 
distribute it in the bulk). Eq. ([1]) becomes {of. Eq. (6) 
in Ref.li): 

w^3cp(B) = w,ep(A) + r {£;,eg(B) - £;seg(Bs)} , (2) 

where A and B denote pure material and material with 
segregant, respectively. Segregation energies are easy to 
obtain in an ab initio supercell approach using the total 
energies of supercells containing grain boundary (or sur- 
face) with and without impurity, £'tot(B) and Etot{^), 
and the same combination for the bulk, ^^(^((B) and 
E^^f{A). For a grain boundary, for example, we have: 



^seg(B) 



gb 



(A)} 

(A)} , 



(3) 



WORK OF SEPARATION AND GHOST 
IMPURITY CYCLE 



The energy release rate (the minimal energy per unit 
area of crack advance) associated with brittle cleavage of 
a grain boundary, Gcicav, is the central quantity charac- 
terising the resistance of the boundary to decohesion in 
the Rice-Thomson-Wang approach. If Gdcav is lower 
than the energy release rate associated with emitting one 
dislocation, C/disi, then cracks remain atomically sharp 
and the crystal breaks in brittle manner. If Gcieav > Sdisi, 
the crack blunts and the crystal is ductile. Impurity seg- 
regation to grain boundaries can either decrease or in- 



where N^^ and N'' denote the number of impurity atoms 
included in grain boundary and bulk supercells, respec- 
tively. Eq. ^ conveys a simple picture in which an im- 
purity atom at the boundary is exchanged with a host 
atom in the bulk; if interstitial positions are involved 
then Eq. ([3]) should by augmented by adding or subtract- 
ing a suitable amount of chemical potentials of the host 
atoms which, again, are the total energies per atom in 
pure bulk. 

Implicit in Eq. ([3]) is that the bulk is sufhcicntly dilute 
so that neither interface nor bulk supercells include any 
additional randomly distributed impurity atoms. In the 
dilute limit, the grain boundary impurity excess F and 
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FIG. 2: Equilibrium structure of the E5{310}[001] symmetric tilt grain boundary in pure copper. Larger and smaller circles 
represent alternating (001) Cu planes. The grain boundary plane contains two inequivalent Cu atoms to which we refer as to 
the "loose" and "tight" sites. The "interstitial" site suitable for segregation of small impurity atoms shown with an asterisk. 



excess volume v'^" per unit area can be found as 



= j{v{B)-Ntno} , (5) 

where V^(B) is the volume of a relaxed grain boundary 
supercell with impurity, A^q is the number of host atoms 
in this supercell and CIq is the atomic volume in pure bulk. 

To separate various aspects of grain boundary weaken- 
ing or strengthening by impurity atoms originating from 
their size, positions, and chemical identity, the "ghost 
impurity cycle" introduced in Ref. |3 and shown in Fig. [l] 
is rather useful. In this cycle, the direct transition 
from unsegregated to segregated state A— >B is replaced 
with a gedanken path through intermediate configura- 
tions A^D— >C— >B for both grain boundary and sur- 
faces, and the respective changes in Wsop are evaluated. 
Configuration C is created from B by removing all impu- 
rity atoms without subsequent relaxation. These missing 
atoms are referred to as "ghosts" since they create forces 
which keep host atoms in place but do not contribute to 
the energy of the system in any other way. Such "ghosts" 
are distorted vacant sites for substitutional impurities or 
centres of expansion for interstitial impurities. Configu- 
ration D is constructed from A in a similar way except 
that we remove the host atoms defined by impurity sites 
in B only if these impurities replace host atoms. For in- 
terstitial impurities, configuration D is not visited (see 
interstitial path in Fig.[T]). 

The same approach is used for the generation of sur- 
face configurations Ag-Dg. Bg and Ag represent the 
equilibrium geometry of free surfaces with and without 
impurity, respectively. Configuration Cg is created from 
Bg by removing impurity atoms, whereas in Dg one re- 
moves only those host atoms that will be replaced with 
impurity in Bg. The remaining atoms in Cg and Dg are 
not allowed to move in response to removal of some of 



their neighbours. Note that the impurity atoms occu- 
pying substitutional positions at a grain boundary can 
be interstitial impurities at free surfaces (or vice versa). 
In terms of Fig. [1] this would mean that the substitu- 
tional path should be used for grain boundaries, whereas 
the interstitial path should be taken for surfaces. More 
generally, one could envisage a situation in which part of 
impurity atoms occupy interstices whereas the other part 
substitute host atoms. We shall describe shortly how to 
deal with such situations. 

As argued in Ref. 8, the change of Wscp at the A^D 
step describes grain boundary weakening due to some 
host-host bonds being broken ("host removal" mecha- 
nism, HR). Transition D-^C corresponds to the distor- 
tion of the atomic structure of pure boundary and surface 
caused by impurity ("substitutional structure" mecha- 
nism, SS). Finally, step C— >B brings in the impurity-host 
chemical interactions and, if relevant, associated changes 
in neighbouring host-host bonds. For oversized impu- 
rity atoms, this step also incorporates the elastic energy 
stored in compressed impurity atoms. These two mecha- 
nisms can not be separated hence we refer to this step as 
"chemical -I- compressed impurity" mechanism (CC). For 
more details regarding the cycle and its implementation, 
the reader is referred to Ref. H. 

The "ghost impurity cycle" outlined above treats in- 
terstitial and substitutional impurities on an equal basis. 
This is vital for the purposes of the present study in which 
we shall be introducing boron into substitutional and in- 
terstitial positions at a grain boundary, sometimes even 
simultaneously. The way to deal with the "co-existence" 
of the substitutional and interstitial paths in Fig. [T] is to 
formally include configuration D into the latter making 
it indistinguishable from A. In such case Eq. ^ can be 
used without making any specific allowances. The same 
applies to surfaces with impurities occupying adatom po- 
sitions. These can be treated similarly to interstitial im- 
purities at grain boundaries. 
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TABLE L Ground-state properties of a— B: equilibrium volume, Vo, rhombohedral lattice constant, ao, rhombohedral angle, 
(f}, internal coordinates of boron atoms, xi,z-i,X2, and Z2, bulk modulus, B, and the energy difference (per atom) between fee 
and a boron, AEfcc, obtained in the present study and other ab initio calculations. Experimental data in the last column are 
from Ref. [ij unless stated otherwise. 



Method 

LDA or GGA? 



FP LMTO PAW 
LDA 



PW-PP LMTO FP LMTO PW-PP PW-PP PW-PP PW-PP US-PP US-PP Exp. 



LDA 



LDA 



LDA LDA LDA GGA GGA 



Reference 


present study 


Ref. 12 


Ref. 13 




Ref. 14 


Ref. 15 Ref. 16 Ref. 17 Ref. 18 


Ref. 11 


Vo (A«) 


6.899 


6.993 


7.05 6.93 


7.06 


6.88 




6.946 7.30 


7.337 [19] 


ao (A) 


4.967 


4.989 








5.034 


4.98 4.967 4.973 




<t> (deg.) 


58.055 


58.063 








58.119 


58.2 58.65 


58.06 


Xl 


0.0106 


0.0105 




0.010 








0.0104 


Zl 


-0.3460 


-0.3458 




-0.344 








-0.3427 


X2 


0.2211 


0.2215 




0.220 








0.2206 


Z2 


-0.3694 


-0.3700 




-0.369 








-0.3677 


B (GPa) 


232 




249 266 


230 


227 




218.4 


224 [19] 


AEf,, (eV) 


1.30 


1.35 


1.43 1.83 


1.31 


1.39 









II. CALCULATION DETAILS 



III. BORON IMPURITY IN BULK COPPER 



The I]5{310}[001] tilt grain boundary is represented 
in our study by a periodic supercell containing two grain 
boundaries with opposite orientation, without any vac- 
uum. Altogether, there are 38 atoms in the supercell, 
with each atom representing one {310} layer, except 
the grain boundary plain which contains two atoms, the 
"tight" and the "loose" sites. These two sites together 
with an "interstitial" site shown in Fig [2l are consid- 
ered as three possible segregation sites for B. Occupation 
of any one, any two, or all three of these sites corre- 
sponds to 0.5, 1, and 1.5 ML coverage in our notation. 
As we shall see in Sec. IIVBI severe atomic relaxation of 
the boundary involving lateral translation of the grains 
may significantly change the local environment of segre- 
gated atoms. We therefore shall be labelling configura- 
tions with respect to positions in which the impurity was 
initially placed. 

To represent free (310) surfaces we use the same su- 
percell with 25 layers of copper, the rest being vacuum. 
Outermost copper layers are replaced with impurity lay- 
ers if we need to model a segregated surface. The situa- 
tion in which a grain boundary containing 0.5 or 1.5 ML 
of impurity cleaves into surfaces with even amount of im- 
purity requires us to double the supercell along the [001] 
direction. Cubic supercells containing up to 108 atoms 
(3x3x3 fee cells) were used to model boron impurities in 
bulk Cu. 

Our first principles calculations employed the full 
potential LMTO method as implemented in the NFP 
code.-SS Calculations were semirelativistic, without spin 
polarisation. We used the local density approximation 
(LDA) in the parameterisation of von Barth and Hedin, 
modified by Moruzzi et alM^ Other parameters (fc- 
point meshes, real space meshes, etc.) were the same 
as in Ref. to which the reader is referred for further 
details. 



A. Pure boron 

Solid boron can exist in a number of relatively sta- 
ble allotropic modifications — rhombohedral, tetragonal, 
and even amorphous. It is not clear which phase corre- 
sponds to the ground state of B at ambient conditions. 
Two rhombohedral phases, a— and /?— B are the most 
likely candidates, a— B becomes unstable at 1200°C and 
converts to f3—B above 1500°C, but (3—B does not trans- 
form back to a— B upon cooling (see [T4| and references 
therein). Hence kinetic effects must impose severe limi- 
tations in this material, a— B was found lower in energy 
by 0.036 eV/atom in Ref. [H and by 0.283 eV/atom in 
Ref. [13, and is assumed to represent the ground state in 
the present study. 

The unit cell of the rhombohedral a— B consists of 
twelve atoms forming an icosahedron. Equilibrium struc- 
ture, bulk modulus and the fee— a energy difference ob- 
tained in our study are in good agreement with other 
calculations and show the usual discrepancies with ex- 
periment associated with the LDA, namely, underesti- 
mation of atomic volumes and consequent overestima- 
tion of the bulk moduli (see Table |T| . Calculations are 
scalar- relativistic, fully relaxed, and employ the 8x8x8 
Monkhorst-Pack mesh of k points. Increasing the k point 
mesh to 12x12x12 changes the total energy by less than 
10~^ Ry, whereas the forces remain within the conver- 
gence criterion lO"'^ Ry/Bohr used throughout the whole 
study. 



B. Copper— boron solid solution 

As noted in the Introduction, the available obser- 
vations do not allow one to conclude unambiguously 
whether the ground state^^ of boron in bulk Cu is inter- 
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TABLE II: Boron in bulk Cu; enthalpy of solution Hs (per 
impurity atom) and the relative dilation volume Jld/^o of in- 
terstitial boron Bi, substitutional boron Bs, and boron dumb- 
ells Bd with different orientation. /S.Hs is the enthalpy relative 
to that of the s(lOO) dumbell (per entity). Nat denotes the 
number of lattice sites in a supercell used in the calculation, 
Q.d is the change of the volume of this supercell when either 
a single impurity (Bi or Bs) or a dumbell Bd is introduced, 
fio is the atomic volume in pure fee Cu, and Ea is the activa- 
tion energy for impurity diffusion. Experimental data on Ea 
are from Ref. 0. Theoretical results are obtained by the FP 
LMTO method unless indicated otherwise. 

Impurity Ngt Hs, eV AH^, eV Q,d/Q.Q Eg, eV 



B, 


32 


1.58 


"0.10 


0.88 


1.31 


Bt 


32 


1.63 




0.93 




Bi 


108 


1.66 


0.06 


0.89 


0.93 


Bs 


32 


1.70 


0.02 


-0.44 




Bs 


108 


1.70 


0.10 


-0.45 




Bd: 












a(lll) 


32 


1.69 


1.70 


0.52 




s(lll) 


32 


1.04 


0.40 


0.38 




a(100)t 


32 










s(lOO) 


32 


0.84 





0.25 




s(lOO) 


108 


0.80 





0.24 





Exp. Bi: 0.57(5) 
B^: 1.15(10) 

^ PAW calculations. 
* Unstable, converts to s(lOO) dumbell during relaxation. 

stitial Bi or substitutional Bg. We calculated the heats 
of solution of the both impurity types using 32 and 108 
atom supercells and find that B^ (in octahedral site) is 
marginally more stable (by 0.04 eV with a 108 atom su- 
percell, Table . The difference is small hence it seems 
reasonable to expect that both B^ and Bs can be found 
at elevated temperatures, as was indeed observed jiSiiS 

However, as shown in Table |TT1 either the insertion 
of a boron atom into an interstice or replacement of a 
host atom at a regular lattice site both lead to signifi- 
cant volume change, fid, negative for Bs and positive for 
Bj. Hence, one could hypothesise that combining Bg and 
Bj into a dimer might eliminate most of the elastic dis- 
tortion of the lattice. To explore this idea, we repeated 
the calculation with B.^ and Bs placed next to each other 
either along the (111) or (100) directions. In addition, 
we also considered two boron atoms arranged symmet- 
rically around a vacant site along same directions. We 
shall refer to the former and the latter as asymmetric 
and symmetric dumbells, respectively. Among these, the 
symmetric (100) dumbell, s(lOO), appears to be the most 
stable (see Table HIl. The heat of solution of the s(lOO) 
dumbell is by a factor of two lower than those of the single 
impurities, indicating that the dumbells should dominate 
at low temperatures and even survive up to the eutectic 
temperature = 1013°C. 

To verify the last point, we use a simple model in which 



boron dumbells B^ together with single impurities B^ and 
Bs are treated as three types of coexisting point defects 
forming an ideal solution. The model is described in 
Appendix A, together with the results obtained within 
this model in the dilute limit. We find in particular, 
that the concentration of dumbells exceeds those of single 
impurities in most conditions unless the temperature is 
close to Te or the boron content is very small. Otherwise, 
all three impurity forms coexist, both in the single phase 
and the two phase regions of the Cu-B phase diagram, 
with the dumbells usually being the dominant kind. 

In fact, this can be anticipated already from the dif- 
ference of the solution enthalpies per entity {Hs for sin- 
gle impurities and 2Hs for dumbells) listed in Table |TT] 
in column AHg. These differences serve to estimate 
the relative amount of defects at terminal solubility (see 
Eq. (jAlSp in the Appendix). 

The fact that adding B to Cu gives rise to three types 
of competing defects is due to a remarkable coincidence 
that B atoms are so perfectly "incompatible" with the 
Cu lattice that the interstitial and substitutional sites 
are nearly degenerate in energy; and furthermore, boron 
dumbells stabilised by this misfit strain turn out to have 
almost the same heat of solution per dumbell as those of 
the single impurities per atom. 

The individual concentrations of the defects might of 
course change if temperature effects, such as atomic vi- 
brations and lattice expansion, are fully taken into ac- 
count. In addition, association of several (three, four, 
etc.) impurities can also play a role, at least at low tem- 
peratures. Nevertheless, we believe that our results pro- 
vide a strong indication that in equilibrium copper-boron 
alloys a significant fraction of the impurities are found in 
"bound" states. 

We are not aware of any metallic system in which dilute 
impurity would aggregate. Interestingly, Dewing noted 
that activity measurements of B in molten Cu suggested 
that B should dimerise in dilute solutions. ■^^ Here we seem 
to arrive at the same conclusion although for different 
reasons as elastic strain does not exist in liquids. Our 
finding might still be relevant to the processing of the 
experimental data on copper-boron melts as the solid 
alloy in these studies is customarily assumed to be an 
ideaP^ or regular-- solution of fully dissociated impurity 
atoms. 

Another observation that stems from the heats of so- 
lution in Table |TT] is the fact that the solubility limits 
indicated in the experimental phase diagram;^ 0.06 at.% 
at room temperature and 0.29 at.% at Tg, are much too 
high in comparison with the enthalpies that we obtained. 
Assuming that impurity atoms are in the form of dimers 
and that boron precipitates as the pure rhombohedral 
a— phase, the above solubilities would translate into the 
Gibbs free energy of solution 0.42 eV/atom at Tg and 0.12 
eV/atom at room temperature. Assuming single impuri- 
ties leads to even larger disagreement. Using the afore- 
mentioned ideal solution model together with the heats of 
solution in Table [Til leads to a three orders of magnitude 
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discrepancy in terminal solubility at Te, which increases 
to more than twenty orders of magnitude at room tem- 
perature (see Appendix A). Temperature effects can be 
noticeable at Tg but cannot explain the discrepancy at 
room temperature. 

Puzzled with this inconsistency, we compared our FP 
LMTO calculations with those by the PAW method^^i^a 
as implemented in the VASP codei^i^^ Hg obtained for 
B; in the 32 atom cell, also shown in Table HIl agrees with 
the LMTO result within 0.05 eV/atom (the difference 
almost entirely comes from the energy difference between 
fee and a— boron obtained by each method, see Table 
Hence we conclude that the heats of solution presented 
in Table [TTl are the correct LDA result. 

The B-Cu phase diagram in Ref. 3 is taken from 
the critical assessment of available experimental data by 
Chakrabarti and Laughlinj^^ in which the data on maxi- 
mal solubility of B in Cu are solely based on experimen- 
tal work by Smiryagin and Kvurtj2i The solubility limits 
mentioned above are those estimated in this latter study 
(0.05 and 0.01 wt.% B translated into at.%) and appear 
to serve more as an upper boundary rather than as exact 
numbers. Chakrabarti and Laughlin indeed comment in 
their assessment that "it is likely that the actual solu- 
bility is even lower than that given by [13] ■" We expect 
it to be significantly lower and appeal to future experi- 
mental work to correct the terminal boron solubility in 
published B-Cu phase diagrams. 



IV. BORON AT A COPPER GRAIN 
BOUNDARY 

A. Work of separation and grain boundary excess 
volume 

The work of separation of a grain boundary at given 
impurity excess is the difference in total energy between 
equivalent pieces of material containing the boundary 
and free surfaces into which the boundary cleaves. Both 
the grain boundary and the surface pieces should be 
taken in their lowest energy state. 

The lowest energy grain boundary can be found by 
comparing the energies of relaxed grain boundaries with 
impurities initially placed into various sites (substitu- 
tional or interstitial). This procedure does not, of course, 
guarantee arrival at the global minimum but is a prac- 
tical alternative to a full optimisation of grain boundary 
structure and includes rigid translation of the grains. 

The surfaces do not require translations, but a com- 
plication here is that one does not know in advance 
the optimal distribution of the impurity atoms between 
two newly created surfaces. Usually the impurity splits 
equally, but not always. More generally, even if equal 
amount of impurity is experimentally detected for a 
cleaved macroscopic sample, there still remains a pos- 
sibility that the surfaces would contain patches with un- 
even impurity coverage. 
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FIG. 3: Convex hull plot for B impurity at Cu(310) surface. 
Circles correspond to the lowest energy surfaces for each cov- 
erage found in our study. 1 ML surface appears to be unstable 
with respect to decomposition into those with 0.5 and 1.5 ML 
coverage (dashed line). 



The combination of surfaces which produces the low- 
est energy at given overall amount of impurity can be 
identified if one employs the convex hull construction. In 
Fig. [3] we plot the total segregation energy to a surface 
G|(,g — rE^^g as a function of impurity excess F. The 
quantity G^^g shows how much the segregation decreases 
the energy of a piece of material containing a surface. 
Therefore, any concave region of the curve indicates that 
there exists a combination of surfaces which provide lower 
energy. In our case, we find that the 1 ML grain bound- 
ary splits into surfaces with coverage 0.5 and 1.5 ML, 
whereas the other boundaries cleave evenly. 

The resulting works of separation are listed in Ta- 
ble mil together with segregation energies and grain 
boundary excess volumes. The maximal Wsep for each 
coverage correspond to the lowest energy grain bound- 
ary and are highlighted in bold. We thus find that at 
0.5 ML boron prefers the interstitial site and at 1 ML 
replacing Cu at both the "loose" and "tight" sites pro- 
vides the best option. Corresponding works of separation 
are compared with those obtained for Bi, Na, and Ag^ 
in Fig. [31 Contrary to these latter impurities, boron in- 
creases Wsep for the whole range of coverages considered. 
The largest strengthening effect is observed at 0.5 ML at 
which Wscp increases by 0.5 J/m^and becomes as large as 
3.81 J/m^. This can be compared with twice the surface 
energy of pure Cu;^ "^Isio = 4.42 J/m^ which provides 
an upper bound for Wsep above which a grain boundary 
would become stronger than bulk. 

Grain boundary segregation energies in Table IIIII com- 
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TABLE III: E5(310)[001] Cu grain boundary with B at various segregation sites: grain boundary excess volume per unit area, 
v^" , average segregation energies per impurity atom, -Escg, and the work of separation, Wsop- Letters "1", "t", and "i" in the 
second line correspond to impurity atoms being initially placed in the "loose" , "tight" , and "interstitial" positions at the grain 
boundary plane (see Fig. [2]) and then relaxed. The optimal impurity distribution between cleaved surfaces is indicated in the 
last line. The quantities corresponding to the lowest energy grain boundaries at each coverage are highlighted in bold. 



Impurity excess, ML 





0.5 


1 


1.5 


2 


Site 




1 t i 


1+t l+i t-l-i 






GB excess volume per unit 
area u^"", A 


0.28 


1.02 0.53 0.45 


0.87 0.78 0.66 


0.98 




Segregation energy Eacg, eV: 
to the (310) surface 
to the grain boundary 






0.40 0.40 0.40 
-0.66 0.40 0.96 


0.17 0.17 0.17 
0.51 0.31 0.29 


0.37 
0.46 


0.30 


Work of separation Wscp, J/m^ 
cleavage mode 


3.35 


2.49 3.35 3.81 
0.5 + 0.5 


3.57 3.24 3.21 
0.25 + 0.75 


3.58 

0.5 + 0.5 





pare favourably with experimental estimations of 0.4-0.5 
eV^ especially for high coverages. These energies assume 
the s(lOO) dumbells to be the ground state of B in bulk 
Cu (see Table |lT| . Had we used interstitial or substitu- 
tional B impurities instead, then the segregation energies 
would have been by 0.8 eV higher. We take this as an 
independent confirmation of the fact that B impurities in 
bulk Cu dimerise. 



Grain boundary excess volumes w^'' in Table IIIII are 
calculated using the dilute bulk limit, Eq. They are 
smaller than those for Bi, Na, or Ag as shown in Fig. [S] 
The dotted hue in Fig.[5]corresponds to the excess volume 
if Cu is notionally considered as an impurity (one could 
think of ^^Cu isotope, for instance). Segregation of such 
"ideal" impurity leaves any grain boundary intact, hence, 
in accord with Eqs. v^''{r) is a straight line with 



r(ML) 

0.5 1 



1.5 




8 10 12 14 



r (nm"^) 



r(ML) 
0.5 1 



1.5 




3 



8 10 12 14 



r (nm"^) 



o 
T3 



FIG. 4: The work of separation of the E5(310)[001] Cu grain 
boundary, Waep as a function of impurity excess F. Data for 
Bi, Na, and Ag are from Ref. The horizontal line corre- 
sponds to twice the surface energy of Cu(310),— 7s ~ 2.21 
J/m^ which provides a natural upper bound for Wsep. 



FIG. 5: Same as Fig. |4] but for the grain boundary excess 
volume per unit area, v^'. The right hand scale gives v^" in 
units of {310} interlayer spacing in bulk copper, do = 0.5576 
A. The dotted line corresponds to a hypothetical "ideal" im- 
purity identical with Cu atoms. 



8 



the slope given by flQ. The fact, that the excess volumes 
of grain boundaries with B lie below the dotted line in 
Fig. [5] indicates a denser packing of atoms at the grain 
boundary with boron compared to the pure boundary. 



B. Atomic structure of the segregated surfaces and 
grain boundaries 

The atomic structure of relaxed grain boundaries with 
0.5 ML of boron is relatively simple. The segregated 
boundaries retain the structure of the equilibrium pure 
boundary shown in Fig. [2] responding to the insertion of 
impurity by either minor shrinking ("tight" site) or ex- 
pansion ("loose" and "interstitial" sites). These bound- 
aries are not shown here. 

The structure of grain boundaries with 1 and 1.5 ML 
are more complex (see Fig. 15]). The exception is the "t-|-i" 
boundary {i.e. boron atoms segregate to the "tight" and 
"interstitial" sites) which does not significantly change 
compared to the pure boundary. However, the average 
segregation energy of boron to the "t-t-i" boundary is 
lower than those for either "t" or "i" 0.5 ML boundaries 
(see Table Uni. 

Among the 1 ML boundaries, the lowest energy corre- 
sponds to the more open "1-l-t" boundary in which boron 
replaces copper at both "loose" and "tight" sites. This 
boundary experiences a large lateral shift of the grains. 
As a result of this shift, the boron atoms lying in adjacent 
(001) planes become nearest neighbours and form boron 
"strings" running along the [001] direction (normal to the 
plane of the drawing in Fig. [S]) . 

The "l-|-i" boundary also entails a rigid translation of 
the grains. This time, however, boron atoms lie in the 
same plane and therefore cannot form [001] "strings". 
The energy of the "1-1- i" boundary is close to the energy of 
the "t-|-i" boundary despite their atomic structures being 
very different. As a matter of fact, the "1-f i" boundary 
can be interpreted as the "i" boundary in which the grain 
boundary plane is shifted normal to itself by one layer 
and another boron atom substitutes a copper atom in 
the adjacent plane. 

The 1.5 ML "l-|-t-|-i" boundary combines the features 
of the "t-|-i" boundary (interstitial boron surrounded by 
six Cu atoms) and the "1-t-t" boundary (boron "strings"). 
The average segregation energy for this boundary is close 
to that of the 1 ML "1-l-t" boundary. 

The following observation is worthwhile here. Substan- 
tial fall off of segregation energy with the amount of seg- 
regant appears to be a common feature of boron doped 
intermetallic compounds. As Lejcek and Fraczkiewicz 
notep. this could be formally described by either in- 
troducing a strong repulsive term into the Fowler- 
Guggenheim segregation isotherm, or by using the stan- 
dard Langmuir-McLean isotherm but with limited num- 
ber of segregation sites. As we observe here, the former 
approach might be physically misleading, at least for the 
Cu-B system. For instance, the 1 ML "1-f t" boundary 



containing neighbouring boron atoms is lower in energy 
than the "t-|-i" boundary where boron atoms are well 
separated. 

The structure of free (310) surfaces with segregated 
boron were obtained by replacing copper with boron in 
the top layer(s) and allowing the surface to relax. For 
fractional coverages 0.5 and 1.5 ML, we tried either to 
substitute half of the host atoms in a layer with impurity 
(substitutional positions) or to place half a monolayer 
of impurity atoms above the top surface layer (adatom 
positions) . The latter resulted in lower energy configura- 
tions. 

During relaxation boron atoms embed into the sub- 
strate, often going beneath the top copper layer. As a 
result, the top copper layer becomes strongly distorted. 
Fig. [3 shows the Cu(310) surface with 0.5 ML of B (and 
0.5 ML of vacancies) in the first layer. The plane of 
B atoms indeed resides below the top layer of Cu atoms 
and just slightly above the second layer, whereas the host 
atoms in the surface layer are strongly displaced towards 
the nearest boron atom. 



C. The reasons behind grain boundary 
strengthening 

We now apply the "ghost impurity cycle" to the grain 
boundaries described in the previous two subsections in 
order to understand why boron segregation has positive 
effect on Wgop- Contributions from HR, SS, and CC 
mechanisms defined in Sec.|T]in terms of the work of sep- 
aration and segregation energies are listed in Tables IIVI 
and|Vl respectively. We choose to evaluate the segrega- 
tion energies in Table |V] assuming interstitial boron B^ 
to be the bulk ground state. With this choice, we can di- 
rectly compare contributions to the SS mechanism with 
those obtained in other studies. Table |Vll The segre- 
gation energies for configurations C and D become just 
energies required to create unrelaxed vacancies at an in- 
terface taken with opposite sign (our convention here and 
in Ref. 8' is that the bulk impurity is always fully relaxed, 
whether this is a real impurity or a vacancy; relaxed "in- 
terstitial vacancy" is just perfect bulk). If one wants to 
change to boron dumbells B^, then the segregation ener- 
gies in Table |V] should be modified as follows, i^sog for 
configuration B decreases by 0.78 eV (the difference be- 
tween the enthalpies of solution of B^ and B^) and i?seg 
for configurations C and D decrease by 1.27/2 eV (half 
the vacancy formation enthalpy in pure bulk, Ref. [sl). 
Works of separation in Table IIVI do not depend on the 
bulk reference as it cancels out in Eq. 

Intuitively, one may expect that segregation of boron 
to interstitial sites at 0.5 ML coverage would reinforce the 
boundary because additional atoms lead to additional co- 
hesion across the interface provided that the boundary is 
not much distorted. Table IIVI supports this expectation 
as the total increase of Wscp by 0.46 J/m^ is provided 
almost exclusively by the CC mechanism. The SS con- 
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FIG. 6: Relaxed E5 copper grain boundary with 1 and 1.5 ML of boron. For 1.5 ML, boron occupies "loose", "tight", and 
"interstitial" sites at the grain boundary plane. For 1 ML, any two of them are occupied (three combinations). Works of 
separation and grain boundary excess volumes corresponding to these boundaries are listed in Table IIIII Larger and smaller 
circles correspond to host and impurity atoms lying in neighbouring (001) planes. 
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FIG. 7: Two side views of the relaxed Cu(310) surface with 0.5 ML of boron. Boron atoms were initially placed instead of the 
top layer of Cu atoms from which every other atom along the [001] direction was removed. During the relaxation, boron atoms 
descend below the next Cu layer causing a noticeable distortion of the latter. 



tribution is also positive but small whereas the HR con- 
tribution for interstitial impurity is zero by definition. 
Two other 0.5 ML configurations, "1" and "t", do not 
increase Wgcp because of the large negative HR contribu- 
tion arising if boron replaces host atoms. 

A similar effect of interstitial boron is found in ab initio 
studies of grain boundaries in Fe, Ni, and Mo (see Ta- 
ble IVip . Boron improves the cohesion at all grain bound- 
aries, and this is mostly due to the CC mechanism. The 



SS contribution enhances grain boundary strength even 
more despite the surface and grain boundary terms be- 
ing themselves negative, i.e. boron distorts free surfaces 
more than grain boundaries. 

As boron segregation proceeds beyond 0.5 ML, the 
boundary is still strengthened although the magnitude 
of the effect is diminished. Higher coverage configura- 
tions necessarily include the removal of host atoms since 
all interstitial sites are already filled at 0.5 ML. Hence, 



10 



TABLE IV: Work of separation Wscp at the vertices of the "ghost impurity cycle" (Fig. [T]) and its change AW^scp due to 
transitions between the vertices. The latter have the meaning of the contributions of the SS, HR, and CC mechanisms. The 
positive sign of AWscp corresponds to cohesion enhancement. 



Impurity Site Excess, Waep, J/m^ AWscp, J/m^ 

ML ~ B C D Total SS HR CC 



i 


0.5 


3.35 


3.81 


3, 


.38 


3 


.35 


0.46 


0.03 





0.43 


1 


0.5 


5) 


2.49 


2, 


.20 


2 


.14 


-0.86 


0.06 


-1.21 


0.29 


t 


0.5 


15 


3.35 


2, 


.87 


1 


.92 


0.00 


0.95 


-1.43 


0.48 


1+t 


LO 


•>7 


3.57 


3, 


.08 


1 


.35 


0.22 


1.73 


-2.00 


0.49 


1+i 


1.0 


5? 


3.24 


2, 


,84 


2, 


.18 


-0.11 


0.66 


-1.17 


0.40 


t+i 


1.0 


?! 


3.21 


3, 


.15 


1 


.96 


-0.14 


1.19 


-1.39 


0.06 


,+t+i 


1.5 


?! 


3.58 


3, 


,33 


1 


.38 


0.23 


1.95 


-1.97 


0.25 


1 


0.5 


?! 


2.15 


1, 


,81 


2 


.20 


-1.20 


-0.39 


-1.15 


0.34 


1+t 


1.0 


?! 


0.96 


0, 


,13 


1 


.38 


-2.39 


-1.25 


-1.97 


0.83 



TABLE V: Same as Table ITVl but for the segregation energies i?sog- The latter are average segregations energies per impurity 
atom (or vacancy). The sign convention is that positive iJscg means that impurity wants to segregate. Segregation energies for 
the pure boundary i5seg(A) are set to zero. The bulk reference used here is an interstitial impurity Bi, not a dumbell Bd as in 
Table Uni 



Impurity 


Site 


Excess, 
ML 










„eV 














ASscg, eV 






B 




C 




D 




Total 


SS 


HR 


CC 


surf 




;b 


surf 


gb 


surf 


gb 


surf 




;b 


surf 


gb 


surf 


gb 


surf gb 


B 


i 


0.5 


1.18 


1, 


.74 


-0.14 - 


-0.10 








1.18 


1 


.74 


-0.14 - 


-0.10 








1.32 1.84 




1 


0.5 


1.18 


0, 


.12 


-0.14 - 


-1.56 


- 


-1.49 


1.18 





.12 


-0.14 - 


-0.07 





-1.49 


1.32 1.68 




t 


0.5 


1.18 


1 


.18 


-0.14 - 


-0.73 


- 


-1.76 


1.18 


1 


.18 


-0.14 


1.03 





-1.76 


1.32 1.91 




1+t 


1.0 


1.16 


1, 


.29 


-0.34 - 


-0.50 


-0.02 - 


-1.25 


1.16 


1 


.29 


-0.32 


0.75 


-0.02 


-1.25 


1.50 1.79 




1+i 


1.0 


1.16 


1, 


.09 


-0.34 - 


-0.65 


-0.02 - 


-0.74 


1.16 


1 


.09 


-0.32 


0.09 


-0.02 


-0.74 


1.50 1.74 




t+i 


1.0 


1.16 


1, 


.07 


-0.34 - 


-0.46 


-0.02 - 


-0.88 


1.16 


1 


.07 


-0.32 


0.42 


-0.02 


-0.88 


1.50 1.53 




1+t+i 


1.5 


1.15 


1, 


.25 


-0.41 - 


-0.42 


-0.03 - 


-0.84 


1.15 


1 


.25 


-0.38 


0.42 


-0.03 


-0.84 


1.56 1.67 


Bi 


1 


0.5 


3.11 


1, 


.63 


1.18 - 


-0.71 


1.20 - 


-0.22 


3.11 


1 


.63 


-0.02 - 


-0.49 


1.20 


-0.22 


1.93 2.34 




1+t 


1.0 


3.01 


1, 


.54 


1.23 - 


-0.75 


1.22 


0.02 


3.01 


1 


.54 


0.01 - 


-0.77 


1.22 


0.02 


1.78 2.29 



the mechanism of strengthening is different. 

To understand this mechanism, it is instructive to com- 
pare boron results with those for bismuth,^ reproduced 
in the Tables for convenience. A striking difference be- 
tween B and Bi (or other oversized impurities studied 
in Ref. H) is the sign of the SS contribution. The SS 
mechanism always increases M^scp for boron (Table IIV[) 
and decreases Wscp for Bi, Na, and even Ag without any 
exception. Compare, for example, boron and bismuth 
at the "loose" site, 0.5 ML coverage. In both cases the 
boundary is weakened, more by Bi, less by B. The neg- 



ative HR contributions are similar, the CC mechanism 
acts so as to strengthen the boundary and is even more 
efficient for Bi than for B. Thus, it is the SS mechanism 
which makes the difference, being positive for B but neg- 
ative for Bi. The "loose" site is not of course the best 
choice for boron, but even here it is much less harmful 
than Bi. The comparison of B and Bi for the 1 ML "1+t" 
case is even more telling. The HR mechanism again, has 
large detrimental effect for both, the CC mechanism acts 
in the opposite direction and is nearly twice as large for Bi 
than for B. Finally, the negative SS contribution is large 
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TABLE VI: Boron impurity at an interstitial grain boundary site in different materials: contributions of tlie SS and CC 
mechanisms in terms of the difference of segregation energies _E|cg — E^^^ (in eV per impurity atom). For the SS mechanism, 
individual surface and grain boundary contributions, Eg^g and E^^g, are also shown. The energies in Refs. [35l[3^l37ii38.39] are 
defined through the difference of binding rather than segregation energies. The results, however, can be compared directly. 
Note that the SS and CC mechanisms in Refs. [s^l and [39[ are referred to as the mechanical and chemical contributions, 
respectively. 



Host Grain boundary 


Method 


Total 


SS CC Cohesion Ref. 






gb - surf surf 


gb gb - surf gb - surf enhancer? 



Fe 


E5(010)[001] 


DMol (LDA) 


1.96 


Fe 


E5(210)[001] 


USPP (LDA) 


0.49 


Fe 


E3(lll)[110] 


FLAPW (LDA) 


1.07 


Ni 


E5(210)[001] 


FLAPW (GGA) 


0.49 


Mo 


E5(310)[001] 


MBPP (LDA) 


2.09 


Cu 


E5(310)[001] 


FP LMTO (LDA) 


0.56 











yes 


[35] 










yes 


[36] 










yes 


[37] 


0.27 


-0.16 


0.11 


0.38 


yes 


[38] 


0.95 


-0.23 


0.71 


1.37 


yes 


[39] 


0.14 


-0.10 


0.04 


0.52 


yes 


present 
study 



enough to make the boundary brittle in the Bi case;^ but 
the large and positive SS contribution for B results in 
the boundary being strengthened. In other words, the 
difference in the distortion pattern of grain boundaries 
and surfaces is itself sufficient to either strengthen the 
boundary or to make it brittle. 

What are the reasons for the SS contribution being 
positive for boron? Let us analyse surface and grain 
boundary contributions for configurations "1" and "1+t" 
for B and Bi in Table |Vl The surface contribution to 
the SS mechanism for Bi is negligible, therefore the neg- 
ative (embrittling) effect comes from the grain boundary 
distortion. For boron, this differs in two ways. Firstly, 
there is always a negative surface term, and secondly, the 
grain boundary term can be large and positive. Even if 
the latter is negative (as in cases "1" and "i"), it is still 
smaller than the surface which renders their difference 
positive. 

It is easy to see why the surface contribution to SS is 
negative for boron. It indicates a sizable distortion of the 
surface region and arises for impurities that can embed 
themselves into surface layers (Sec. lIVBp . This would 
be possible for small impurities, especially if they prefer 
interstitial positions in the bulk. 

It is less obvious why the grain boundary contribu- 
tion to SS tends to be positive. SS contribution in the 
"ghost impurity cycle" is the energy change when a pure 
grain boundary with preinserted unrelaxed vacancies D 
is further deformed as prescribed by "impurity ghosts" 
to arrive at configuration C. Atoms in configuration D 
would want to relax towards the vacancies, whereas de- 
formation corresponding to large impurity atoms forces 
them to move further away. The total energy increases 
and the SS contribution is negative (embrittling). For 
small substitutional impurities this is reversed — during 
the D— >C transition atoms move towards the vacancies. 
Hence, the energy decreases and SS is positive (cohesion 



enhancing). If boron occupies both interstitial and sub- 
stitutional sites, atomic displacements are more complex. 
However, the fact that the grain boundary excess volume 
is always smaller than that of the "ideal" Cu-like impu- 
rity (Fig. [5]) indicates that on average the grain boundary 
shrinks rather than expands. 

The above reasoning relies only on the property of 
boron atoms being "smaller" than host atoms and there- 
fore seems applicable to other undersize impurities, at 
least for the light metalloid impurities. Boron was found 
to reinforce grain boundaries in all materials studied (Ta- 
ble IVI[) . Carbon segregation increase o'^^i"^° or slightly 
decreases VFsopr^ whereas H, N, and O weaken grain 
boundaries .i^^^i^ In the latter case, the embrittling 
propensity is due to the CC contribution, which becomes 
large and negative. Janisch and Elsasser— suggest that 
this should be the case for light species whose outer elec- 
tronic shell (Is for H, 2p for N and O) falls near or below 
the bottom of the valence band of the host metal. 

Negative CC indicates that the insertion of an impu- 
rity atom into a prepared hole (configuration C) weakens 
atomic bonds across the interface. This could be the case 
if the impurity affects the bonds between the neighbour- 
ing host atoms by means of withdrawing electronic charge 
from them — which is known as the electronic mechanism 
of embrittlement. An alternative explanation, advocated 
in Ref. [H, is Cottrell's "unified theory" which refers to 
the position of the impurity levels relative to the Fermi 
energy of the host metal42. According to this theory, 
interstitial impurities whose valence electrons lie close 
to the Fermi level, would form predominantly covalent 
bonds with host atoms and hence prefer grain boundaries 
over surfaces due to a higher coordination in the former 
(the Cottrell ^/z factor). That means positive CC and 
cohesion enhancement. On the other hand, impurities 
with valence states lying high above or deeply below the 
Fermi level would form polar bonds with the host atoms 
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and turn into screened ions, for which the surface envi- 
ronment is more favourable. This results in a negative 
CC contribution which weakens the boundary. 

In our recent calculations of Cu grain boundary with 
inert gas atom impurities (He and Kr) we also observed 
a large negative CC contribution leading to catastrophic 
embrittlement4S As no charge transfer to or from inert 
gas atom is expected, the embrittling effect in this case 
must be related to the Pauli exclusion principle. One 
may, therefore, hypothesise that a similar mechanism can 
act for impurities with nearly completed p-shell, such as 
fluorine, oxygen and, to lesser extent, nitrogen, as filling 
the impurity shells with metal electrons would effectively 
render the dopant atoms inert gas like. The question as 
to how significant this "inert gas atom" mechanism is in 
comparison with others requires a separate investigation 
and is outside the scope of the present study. 



V. CONCLUSIONS 

The effect of boron impurities at the S5(310) [001] grain 
boundary, (310) surface and in the bulk of Cu is inves- 
tigated on the basis of first principles calculations using 
the full potential LMTO method. 

1. We find that B strengthens the boundary in the 
whole range of coverages studied (up to 1.5 ML) with 
the maximal effect achieved at 0.5 ML. Combined with 
the observed ability of B to remove harmful impurities 
such as Sb from the copper boundary,— this makes boron 
a particularly attractive alloying addition. 

2. The reasons behind grain boundary strengthening 
at 0.5 ML and higher coverages are different. 0.5 ML 
corresponds to all interstitial positions at the boundary 
being filled by boron atoms providing therefore addi- 
tional cohesion between the grains while not distorting 
the boundary much (the CC mechanism). 

3. At 1 and 1.5 ML boron begins to substitute host 
atoms at the boundary leading to significant distortions 
and lateral translations of the grains. The SS contribu- 
tion, however, remains positive and acts so as to increase 
Wscp- We demonstrate that the difference in the sign of 
the SS contribution proves to be solely responsible for 
the opposite effect of B compared to embrittling species 
such as Bi. 

4. Distortion of a free surface by segregated boron 
atoms further increases Wscp, but it is not a decisive fac- 
tor. 

5. Introducing boron into bulk Cu leads to a peculiar 
situation in which substitutional and interstitial impuri- 
ties are rather close in energy. Combined together, they 
form a strongly bound dimer held by elastic forces of 
the host lattice. Remarkably, the heat of solution of the 
lowest energy s(lOO) dumbell (per dimer) is also close to 
the heat of solution of boron single impurities (per atom). 
Thus a sizable proportion of boron atoms should be found 
in a bound state in most experimental conditions, even 
at high temperature. 



6. A large discrepancy between calculated heats of 
solution and experimental estimations for terminal solu- 
bility of B in Cu is discovered. We are inclined to think 
that the solubility limits suggested in Ref. 34 and then 
translated into existing B-Cu phase diagram, are over- 
estimated by a few orders of magnitude and hope that 
our findings inspire experimental work on the updated 
version of the phase diagram. 
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APPENDIX A: EQUILIBRIUM 
CONCENTRATION OF COEXISTING 
IMPURITY TYPES IN AN IDEAL SOLUTION 

We describe a thermodynamic approach which we use 
in the paper to estimate the equilibrium concentration 
of boron single impurities (B^ and B^) and dimers (B^) 
in bulk Cu. The approach follows the canonical treat- 
ment of ideal solid solutions proposed in Ref. 44 and is 
applicable to any binary system in which the heats of 
solution of the impurity species occupying an interstitial 
position, a substitutional position, or forming a dumbell, 
are comparable. 

The reader may be surprised by the fact that the 
model outlined below ignores thermal vacancies alto- 
gether. Indeed, at a first glance this looks inconsistent 
given that the heats of solution of boron in Cu listed in 
Table |lT] are comparable with the vacancy formation en- 
thalpy HJ =1.27 eV.— We omit vacancies deliberately as 
the equilibrium impurity concentrations do not depend 
on the vacancy formation enthalpy, hence they do not 
change even if no vacancies at all are allowed. Indeed, if 
a crystal contains thermal vacancies with concentration 
c^, an additional term will appear in Eqs. (jA4p . (|A5p . 
and (jA7p . and, consequently, one more relation will be 
added to the system of equations (|A8|) . However, c„ can 
be eliminated from (jA8[l explicitly leading to exactly the 
same set of equations (jA9[) - (jAll[) as below. The equilib- 
rium impurity concentrations, obtained as the solution 
of Eqs. (|A9p - (|All[) . will therefore not depend on the va- 
cancy concentration either. (Note in passing that the 
reverse is not true, i.e. the equilibrium concentration of 
vacancies does depend on impurity concentrations.) 
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1. Variables and definitions 

Consider a large piece of crystal Ai^^Bx containing 
N lattice sites with tt-q host atoms A and nt, impurity 
atoms B. The latter in turn, include ris substitutional 
impurities, interstitial impurities, and Ud dumbells: 

ns + Ui + 2nd = nb ■ 

If the crystal is sufhciently large so that any surface ef- 
fects can be neglected, the resulting equilibrium concen- 
trations should depend on and rif, only through the 
composition of the solid solution 



rib 



(Al) 



Concentrations of host atoms (ca) and impurity of any 
type (cs,Ci, and c^), defined "per lattice site" here, must 
satisfy the following two constraints: 



NCa = Ua 

N {cs + Q + 2q) = Ub , 



(A2) 
(A3) 



indicating that the total number of atoms of each species 
is conserved. In addition, we have the "site balance" 
condition as every lattice position should be occupied by 
either a host atom, an impurity atom, or an impurity 
dumb ell: 



Q = 1 . 



(A4) 



Overall, there are five variables (A^, Ca, Cs, q, and Cd) and 
three constraints (|A2p - (|A4p . hence the system has two 
degrees of freedom. 

If defects do not interact, the total energy of the system 
is linear in defect concentrations: 



E ^ N (CaEa + Cs£s + C.,£i + 2Crfe<i) 



(A5) 



where Eq is the energy (per atom) of the pure crystal, 
whereas £s,£i, and £d are the energies per impurity atom 
defined by Eq. (|A5|) . In practice these are usually found 
from total energy calculation of supercells containing a 
single defect of each type. Enthalpies of solution iJ^, 
such as those listed in Table [III are related to the e's in 
a simple way: 



Ht — Si 



Hg — £s — Hs — £b : 



where Sb denotes the energy per atom of species B in its 
pure state. 

At zero pressure, the Gibbs free energy of the crystal 

is 



G = E-TS, 



(A6) 



where T is the temperature and S is the (configurational) 
entropy given by 



S 



-kN 



Ca log Ca + Cs log Cs + Cd log Cd - Cd log T] 



+ Ci logCi -I- (a - Ci) log (a - Ci) - aloga , (A7) 



where k is the Boltzmann constant, rj is the number of 
equivalent orientations of the dumbell, and a is the num- 
ber of interstitial sites per one lattice site. For the octahe- 
dral interstices in the fee lattice a = 1, whereas 77 = 3 for 
the (100) dumbeU in cubic crystals. Eqs. ((X5|) and (K7\i 
do not take into account atomic vibrations, but these 
can be easily included, for example, at the level of quasi 
harmonic approximation!^ 



2. Equilibrium concentrations 

The equilibrium impurity concentrations are those that 
minimise the Gibbs free energy (|A6[) subject to con- 
straints (|A2|) - (|A4p . The minimisation leads to the fol- 
lowing system of five equations 

G - Hana - fJ-bnb = 

e,-l-fcr(l + logCa)-/ia- A/A ^ 

es + kT{l + logCs)- fib- X/N = (A8) 

2£d + fcT(l + logCrf/ry)-2Aib-A/A = 

Si + kTlog[ct/{a - c^)] - fib = 0, 

where fia and fib are Lagrange multipliers associated with 
Eqs. (|A2[) and (|A3[) and therefore have the meaning of 
the chemical potentials of species A and B, respectively. 
A is the Lagrange coefficient corresponding to Eq. (jA4p . 
Eqs. together with constraints (|X2|) - ((X4)) are suffi- 
cient to determine all the unknowns. 

Elimination of fia, fib, and A leaves the following two 
independent relations containing only the impurity con- 
centrations: 



Ci/a 



(1 - c,/a)"+ic, 
{ci/a)cs 
(1 - Ci/a)cd/rf 



exp 
exp 



(£^ " £s) 

kT 

{Sj +£s'- 2£d) 

kT 



(A9) 
(AlO) 



in which the reader might immediately recognise "quasi 
chemical" relations describing defect reactions. Eq. (jA9p . 
in particular, corresponds to the conversion of a sub- 
stitutional impurity into an interstitial one, whereas 
Eq. (jAlOp describes the dissociation of a dumbell into 
an interstitial and a substitutional impurity. 
Eqs. ((X9|l - (|A10|) together with the relation 



-f- (1 — x)ci + (2 — x)cd = X 



(All) 



can be used to find all three concentrations Cs,Ci, and 
Cd- [Eq. (jAlip readily follows from constraints (|M|-(|^4| 
and the definition of x (jAip j. 

The above consideration applies to an ideal solid so- 
lution of arbitrary composition x. If the alloy is dilute, 
x <C 1, finding the solution of Eqs. (|A9p - (jAlip simplifies 
and reduces to solving the quadratic 



[27f 



C-di 



a{a + 2)1 + {csi + a + 2x)z- 



. = 0, 
(A12) 
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where Cdi and Csi denote the following exponentials 

Tjd Tji 



exp 



exp 



kT 



kT 



exp 



= exp 



kT 



kT 



Equilibrium concentrations are then obtained using the 
positive root of (|A12p . as I 



_ 2 2 

Cd — V e^i ; 



z+ (A13) 



3. Solubility limit 

As pointed out in [3] , the advantage of the above ap- 
proach is that it produces not only the equilibrium con- 
centrations of defects but also the chemical potentials of 
the species. The chemical potential of the solute, in par- 
ticular, can be restored from the equilibrium concentra- 
tions Cd, Ci, and Cs using any of the following equations: 



fib = es + fcT [logCs -I- alog(l - c,;/a)] 
= e^ + kT\og[ci/{a-c,)] 

= Ed + ■^kT [log Cd/v + alog{l ~ Ci/ a)] 



(A14) 



These three relations simply express the fact that im- 
purity atoms participating in any of the three types of 
defects considered here, namely interstitial impurity, sub- 
stitutional impurity, and impurity dumbells, are in equi- 
librium with each other, hence their chemical potentials 
must be equal. 

Once the chemical potentials are known, it is straight- 
forward to find the limiting solubilities by considering 
the equilibrium between the A-rich and the i3-rich phases 
with terminal compositions. The terminal compositions 
are those that make /i^ and fib in the both phases equal. 

For our purposes, however, it is sufficient to assume 
that the i3-rich phase is a pure B crystal, such as the 
rhombohedral a-boron, with Hb — Sb- Using again the 
dilute limit, from (|A14[) we obtain the maximal defect 
concentrations in Ai^xBx as 



Cs = exp 

= a exp 
Cd = Jyexp 



El 
kT 

kT 
_2Hl 
kT 



(A15) 



These, according to (jAlip . define the limiting solubility 
X™ as 



(A16) 
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FIG. 8: Equilibrium concentration of Bd, Bi and Bs in 
Cui-^B^ at T = 1000 K given by Eqs. (|Il2)l - (|Il3)l as a func- 
tion of boron content x. Tfie vertical dotted line corresponds 
to tlie solubility limit estimated according to Eq. (|A16p . 



4. Results for boron in copper 

Fig. [8] shows the equilibrium concentrations of B;, B^, 
and Bd at T — lOOOK as a function of the boron content 
X. These were obtained using Eqs. (|A12p and (jA13p with 
a = 1,77 = 3, and the enthalpies of solution Hg listed 
in Table [Til (those calculated in a 108 atom supercell). 
The limiting solubility x™ given by Eqs. (|A15p - (|A16p is 
shown with a vertical dotted line. 

According to Fig. [51 the dumbells strongly prevail at 
and above a;™. (Supersaturated solutions may arise if 
the excess boron precipitates as some metastable phase, 
higher in energy than the rhombohedral a-B. Such a 
scenario, however, is not supported by experimental 
observations-^'). With x decreasing, the loss of config- 
urational entropy should eventually overweigh the ener- 
getic advantage of forming dumbells, giving rise to the 
crossover between concentration of dumbells and single 
impurities. This is indeed observed in Fig. [SI although 
the crossover concentrations seem too small to be exper- 
imentally detectable. 

Fig.[51shows the temperature dependence of the impu- 
rity concentrations at fixed x in the form of the Arrhenius 
plot logc = f{l/T). Concentration x for this plot is cho- 
sen as the limiting solubility in Fig. [51 As a result, the 
curves in Fig. [51 have a kink at T = 1000 K which cor- 
responds to the precipitation of the second phase. The 
curves to the left of the kink are the equilibrium defect 
concentrations in a single phase crystal Cui-^^B^; given 
by (|A12p - (|A13p . whereas concentrations to the right of 
the kink correspond to a two-phase equilibrium between 
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0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 
1000/T(K'^) 

FIG. 9: Equilibrium concentration of Bd, Bi and Ba and 
the total amount of boron, x = Ci + Cs + 2cd, as a function 
of (the Arrhenius plot). The alloy composition is taken 
as the terminal solution at T = 1000 K (vertical dotted line 
in Fig. |8}. The kink on the curves corresponds to the pre- 
cipitation of the second phase (assumed to be pure a— B). 
The vertical dotted line indicates the eutectic temperature 
Te = 1013°C. 



Cui-x^x and pure a-boron. The latter are given by 



Eq. (|A15p and are just straight lines in the Arrhenius 
coordinates. 

We again observe the dominance of boron dumbells in 
the whole range of temperatures, except the narrow re- 
gion close to the eutectic temperature Tg = 1013°C (ver- 
tical dotted line in Fig. [9]). The crossover point appears 
here for the same reason as in Figs. [5] and rapidly moves 
to higher temperatures with x increasing. If, for exam- 
ple, for this plot we use x corresponding to the limiting 
solubility at T^, then the crossover would move above 
both the eutectic temperature and melting temperature 
of pure copper — 1085°C, so the dumbells would 
dominate everywhere. 

To summarise, we observe that boron dumbells are 
the essential component of Cui_2;Ba; solid solutions both 
in the single phase and two phase regions of the Cu-B 
phase diagram. As a matter of fact, the concentration of 
boron dumbells B^ exceeds those of single impurities B^ 
and Bg in most conditions. High temperature and small 
boron content tend to make these concentrations compa- 
rable at best. The dumbells can be suppressed only in 
very diluted samples where the impurity concentrations 
are likely to fall below the detection limit anyway. 

The terminal solubility of B in Cu appears to be lower 
than that indicated in published phase diagramsi^ It is 
of the order of a few ppm at Te (cf. 0.29 at.% in 0]), and 
is as low as 10~^^ at room temperature (cf. 0.06 at.% 
in 0). The fact that the second phase precipitates as 
a solid solution of Cu in B^i22 rather than as the pure 
boron (assumed here) should lead to even lower limiting 
solubilities. 
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